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Abstract. We calculate the shear viscosity r\ in the quark-gluon plasma (QGP) phase within a virial 
expansion approach with particular interest in the ratio of rj to the entropy density s, i.e. rj/s. The virial 
expansion approach allows us to include the interactions between the partons in the deconfined phase and 
to evaluate the corrections to a single-particle partition function. In the latter approach we start with 
an effective interaction with parameters fixed to reproduce thermodynamical quantities of QCD such as 
energy and/or entropy density. We also directly extract the effective coupling ay for the determination of 
77. Our numerical results give a ratio rj/s ~ 0.097 at the critical temperature T c , which is very close to the 
theoretical bound of 1/(4tt). Furthermore, for temperatures T < I.8T C the ratio rj/s is in the range of the 
present experimental estimates 0.1 — 0.3 at RHIC. When combining our results for rj/s in the deconfined 
phase with those from chiral perturbation theory or the resonance gas model in the confined phase we 
observe a pronounced minimum of rj/s close to the critical temperature T c . 

PACS. I2.38.Mh Quark-gluon plasma - 25.75.Nq Phase transition in Quark-gluon plasma - 21.65.Qr 
Quark matter/nuclear matter - 51 .20.+d Viscosity, diffusion, and thermal conductivity 



1 Introduction 

The exploration of the phase structure of QCD is a chal- 
lenging task for modern theoretical physics. Lattice QCD 
(1QCD) calculations for vanishing quark chemical poten- 
tial fj, q have shown chiral and deconfinement phase tran- 
sitions at a critical temperature T c of about 150 to 200 
MeV Q3Q]- New methods to extend 1QCD calculations to 
finite chemical potentials fj, q have been developed in the 
past decade, i.e. a multi-parameter reweighting [HUMS], a 
Taylor expansion around /i q ~ [HEIE] and an imaginary 
chemical potential method [51ITU]; however, their validity 
is limited to a region /x q < T [11] . At finite net quark den- 
sities presently some modeling of QCD is needed in order 
to explore the QCD phase diagram at least on a quali- 
tative level. Indeed, effective approaches indicate a very 
rich phase structure of QCD, i.e. 1) a quark-gluon plasma 
(QGP), 2) a color superconducting phase, 3) a color-flavor 
locking phase and different further combinations (cf. Refs. 
[T2llT3lH^ll5U16IITT] ). The explicit phase structure, how- 
ever, is model dependent (except for very high /i q ) and 
not accessible by experiment in the full (T, /z q ) -plane. 
Experimental information can be qualitatively obtained 
from astronomical observations and more quantitatively 
from relativistic heavy-ion collisions. In the last decade 
such experiments have been performed at the Relativistic 
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Heavy-Ion Collider (RHIC) at Brookhaven National Lab- 
oratory (BNL) and even long before - at lower energies - 
at the Super-Proton-Synchrotron (SPS) at CERN. 

Results from 1QCD calculations have early been inter- 
preted as indicating a weakly interacting system of par- 
tons since the entropy density s and the energy density e 
are close to the Stefan Boltzmann (SB) limit for a rela- 
tivistic noninteracting system, i.e. only « 10 — 15% lower. 
However, more recent experimental observations at RHIC 
have drastically changed the notion about the QGP. In 
this context, one of the most intriguing experimental find- 
ings is the large elliptic flow of hadrons at RHIC [18, 19 
which is significantly larger than at SPS energies [2"UI|21| . 
In particular the transverse momentum dependences 
of the elliptic flow «2(pt) at RHIC is close to predictions 
from non-dissipative hydrodynamical simulations 22,23, 
HM1I25II26] around midrapidity (|y| < 1). This result has 
led to the BNL announcement about the discovery of the 
nearly perfect fluidity of the strongly-coupled quark-gluon 
plasma (sQGP) [27 p g j produced at RHIC. Especially, due 
to the asymptotic freedom of QCD, the prevailing idea has 
been an expectation of large shear viscosities in a weakly 
interacting QGP (wQGP) at very high densities and/or 
temperatures. Because of the evident failure of these as- 
sumptions at RHIC conditions the novel notion of a strong 
quark-gluon plasma (sQGP) has been put forward to char- 
acterize the observed strong coupling properties of the 
QGP close to (or slightly above) the critical temperature 
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T c that keep viscous effects low at RHIC. Therefore, the 
inclusion of interactions between partons is mandatory to 
consistently describe the QGP in the region close to T c . 

In going beyond the standard dynamical quasiparti- 
cle picture [16,29,301 1311132"] - incorporating interactions 
in terms of a width in the spectral functions of the par- 
tons - we have recently developed a generalization of the 
classical virial expansion approach to calculate the QCD 
partition function in the partonic phase with an interac- 
tion inspired by lattice calculations [33] . We have obtained 
an Equation-of-State (EoS) for the partonic QGP that is 
well in line with recent three-flavor QCD lattice data [T] 
for the pressure, speed of sound and interaction measure 
at nonzero temperature and vanishing chemical potential 
(fi q = 0). Since in the virial expansion approach all ther- 
modynamic quantities are based on an explicit parton in- 
teraction in form of a potential, this approach is also the 
ideal starting point for a consistent description of dynami- 
cal properties of the QGP like the shear viscosity rj, which 
is the stationary limit of a nontrivial correlator [54ll35| . 

In this work we focus on the calculation of the ra- 
tio r)/s, well known as specific viscosity, within the virial 
expansion approach. It was shown some years ago for su- 
persymmetric Yang-Mills (SYM) gauge theory using the 
Anti de-Sitter Space/ Conformal Field Theory (AdS/CFT) 
duality conjecture [3S] that 
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which is denoted as KSS bound. This limit is close to the 
simple bound rj/s > 1/15 from the kinetic theory uncer- 
tainty principle [37] ■ Additionally, it has even been spec- 
ulated that this bound might hold for all substances [55] . 
Furthermore, it has been found for atomic and molecular 
substances that the specific viscosity exhibits a minimum 
in the vicinity of the liquid-gas critical point; this also sug- 
gests the possibility of a minimum of rj/s for QCD at the 
critical temperature T c which has been already noted in 

Refs. [3H1E3 



2 Kinetic theory 

We start by recalling results well known from literature. 
Here kinetic theory is a convenient framework to start the 
investigation of the viscosity coefficient rj [35]. In an ul- 
trarelativistic quark-gluon plasma, i.e. where the temper- 
ature T is much larger than the constituent masses rrii, 
the estimate for the ratio of the shear viscosity 77 and the 
entropy density s is [3 7) 



s s 15 



(2) 



where pi is the local density of quanta i transporting on 
average a momentum (p)i over a momentum-degradation 
length (mean- free-path) A^. Following the more detailed 
kinetic theory derivations of Ref. [40] , the familiar factor 
I of elementary non-relativistic kinetic theory has to be 
replaced in the ultrarelativistic domain by 4/15. 



Evidently, two aspects play a crucial role in the expres- 
sion ([2]) for rj/ s. The first one is the choice of the equation- 
of-state for the determination of the entropy density s. In 
general - for such calculations - an ideal gas EoS is used. 
This assumption might be justified for high temperatures 
T ;§> T c , where the deviation of the QGP equation-of-state 
from the ideal gas limit might be eventually neglected. On 
the other hand, close to the critical temperature T c , this 
approximation should be abandoned due to large devia- 
tions from the SB limit. 

The second crucial point is given by partonic dynam- 
ics itself: The effects of the interaction between the con- 
stituents directly determine the mean-free-path A and its 
calculation is the aim of several investigations within dif- 
ferent approaches. Following Ref. [H] the mean- free-path 
is the inverse of the interaction rate for nearly massless 
quanta J 1 ,, i.e. A; = 1/-Ti, which can be calculated to low- 
est order in the coupling constant g from the imaginary 
part of the quark and gluon selfenergy. Another possibil- 
ity is to relate the mean-free-path to the transport cross 
section a t by A.; = (pert)" 1 . This quantity - in a dense 
partonic system - may be related to the transport pa- 
rameter q governing multiple scattering |421I43| or can be 
directly modeled in a Debye-screened form [44] ; the latter 
leads to an analytic expression for at [55] . Using this last 
formulation and additionally adopting the approximation 
(p) ~ 3T the specific viscosity becomes 
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In particular, assuming that the elastic gluon scattering 
matrix element in a dense partonic medium can be mod- 
eled by a Debye screened interaction [4~4l[4~5"] the relevant 
transport cross section reads 



<r t (f) = J da, 



1 sin 



= a 4z(l + z) [{2z + 1) ln(l + l/z) - 2] , (4) 

with the total cross section ao(s) = 9ira v (s)/2jil CI . Here 
ay = ay(T) and /i scr are the effective temperature de- 
pendent coupling constant and the screening mass, respec- 
tively, and z = /i 2 cr /s. For simplicity, we will assume ao 
to be energy independent and neglect its weak logarith- 
mic dependence on s in the relevant energy range and set 
s fts 17T 2 . We recall that the sin 2 (9 cm weight arises in the 
transport cross section because large-angle scatterings are 
most effective in momentum degradation [37j . The cross 
section Q is a monotonic function of /i scr , which plays a 
crucial role for the results of a t . It is important to em- 
phasize that the coupling constant and the Debye mass 
are not independent parameters in the calculation because 
fi scr is determined by the value of the coupling ay itself. 
In perturbation theory one explicitly obtains 
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(5) 



in gluon dynamics. 

For numerical estimates a specific form for ay(T) has 
to be chosen. As a quantitative reference the long-range 
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part of the strong coupling constant extracted from the 
free energy of a quark-antiquark pair in 1QCD |46j has 
been used in the literature [57] ■ We point out that any 
extraction of a coupling constant ay(T) from 1QCD is 
model dependent and deviations (or agreement) of any 
parametrization from (with) the lattice data have to be 
considered with care. Accordingly, several parametriza- 
tions can be found in the literature [29,16,32,46,47.. 

At this point, we emphasize that in previous works [47] 
these two aspects - equation of state and transport cross 
section - are uncorrelated. Therefore a consistent approach 
is needed which provides i) a realistic equation-of-state - 
and thus a proper entropy density - as well as ii) a trans- 
port cross section within the same framework. 



3 The virial expansion approach to QCD 

As pointed out in the previous Section, a consistent ap- 
proach requires two fundamental ingredients: i) a calcu- 
lation of the equation-of-state including the interactions 
between the constituents and ii) an extraction of the ef- 
fective coupling that enters the estimate for the trans- 
port cross section. Both requirements can be achieved 
within the virial expansion formalism developed in a pre- 
vious work [33] . where a detailed derivation of the parti- 
tion function Z(T,V), of all thermodynamic quantities - 
such as pressure, entropy density, interaction measure and 
sound velocity- and of the EoS of the QGP at vanishing 
and finite ji q has been presented. We achieve an expansion 
of In Z in powers of the logaritm of the partition function 
in the Stefan Boltzmann limit lnZ^ ), i.e. (cf. [55] ) 



lnZ«lnZ(°) + |(lnZW 
with the second virial coefficient 

b 2 = J d 3 r (e-^« - 1 

The pressure is simply obtained as 
P = Tin Z. 



(6) 



(7) 



(8) 



The other quantities can be calculated from the pressure 
P or the partition function Z using thermodynamic rela- 
tions. In particular, for the entropy density one obtains 

M 

dP 

s = W 0) 

For an application of this formalism to the QGP a specific 
choice of the interaction potential W\2 has to be made. 
Following Ref. [33] we use an effective quark-quark poten- 
tial inspired by a phenomenological model which includes 
non-perturbative effects from dimension two gluon con- 
densates (that reproduce the free energy of quenched QCD 
very well) . The effective potential between the quarks ex- 
plicitly reads 



W ia (r,T)= -- 



7T 1 

12 r 



C-2 



2N r T 



-M(T)r 



(10) 



where C2 is the non-perturbative dimension two conden- 
sate and M(T) a Debye mass estimated as 



M(T) = V^c/3 + JVf/6 gT = gT, 



(11) 



where we have neglected any scale dependence in the cou- 
pling constant. By using the potential fj 10[) we assume the 
interaction to be the same for quarks and antiquarks and 
neglect explicit gluon contributions. The latter are en- 
coded in parametric form in the interaction (|10|) . In short, 
we generalize the Yukawa-liquid model for the QGP in- 
vestigated in Refs. [451145] before. For a detailed explana- 
tion of this interaction we refer the reader to Ref. [33] . A 
comparison with three-flavor 1QCD calculations with al- 
most physical masses from Ref. 1 shows that a coupling 
parameter g — 1.30 allows for a good description of all 
thermodynamic quantities in the temperature range from 
0.8 T c to 5 T c . A detailed discussion of the validity of the 
virial expansion truncated at the second term is given in 
Ref. [33]. 

In order to calculate a transport cross section with this 
interaction the coupling ay has to be extracted from V\ . 
Following Ref. [5U] we define the coupling in the so-called 
qq-scheme, 



7t dr 



(12) 



The coupling a qq (r, T) then exhibits a maximum for fixed 
temperature at a certain distance denoted by r max . By 
analyzing the size of the maximum at r max we fix the 
temperature dependent coupling, ay(T), as 



a v{T) = a qq (''max,7 1 ) . 



(13) 



Before we start with the calculation of rj/s two important 
aspects have to be pointed out: first we only consider the 
contribution of our dynamical degrees of freedom (quarks 
and antiquarks) for rj/s, since here the gluons are mass- 
less and interaction free with respect to each other. The 
fermion-gluon interaction is only space-like and included 
in the potential flTTj]) . Therefore, the gluons do not con- 
tribute to the ratio rj/s, in contrast to QCD perturbation 
theory, where the contribution of the quarks and of the 
gluons are roughly of the same order [5T1I52] . We thus 
consider the quark specific viscosity, i.e. 



rj AT 



(14) 



s q 5 s q a t 

where the quark contribution of the entropy density s q is 

s q = s - Sg, (15) 

with the gluon contribution to the entropy density in the 
SB limit given by 



32tt 2 



(16) 



The second important aspect is the applicability of Eq. (|4]) 
to calculate the transport cross section. We recall that this 
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expression for ot has been derived for a Debye screened in- 
teraction [44], whereas the effective potential W12 contains 
not only such a term but also a purely screened part ~ C2 ■ 
A preliminary analysis, which actually has lead to the final 
formulation of Ref. |33j , shows that the main contribution 
to the virial coefficient is given by the Coulomb screened 
part of the potential. This confirms the observations of 
Ref. [53], that the exact value of C2, which modulates the 
purely screened part of the potential, is not so important 
for a good reproduction of the free energy. Accordingly, we 
may use Eq. ([?]) as a good approximation for the transport 
cross section. 

The advantage of this formalism is that for the Debye 
screening mass no further approximations must be done, 
as, for example, the assumtion of the validity of the pertur- 
bation theory that leads to the specific expression given 
in Eq. §5§ and is used often in the literature [4"Tll48j . In 
our model, the mass is directly given by Ea. (|lll) and in 
this way is calculated independently of the coupling ay 
but, at the same time, consistently -in a thermodynamic 
sense- with the properties of the QGP because the param- 
eter g = 1.3 is fixed by thermodynamic quantities. This 
demonstrates the important difference to other effective 
approaches used to understand the enhancement of the 
parton transport cross section, that was concluded from 
the elliptic flow measurements at RHIC. 

We, furthermore, point out an analogy between strong 
electromagnetic and strong quark-gluon plasma coupling 
presented first in Ref. [IS]: The Debye screening length, 
i.e. ro = cannot be used as a cut-off parameter 

for the calculation of at because the interaction range is 
larger than ru- Studies for complex non-relativistic elec- 
tromagnetic plasmas show the importance of these effects 
as demonstrated in Ref. [54] . Furthermore, by assuming 
a non-screened Coulomb interaction a modified Coulomb 
logarithm A* was derived, which leads to a sizeable en- 
hancement of the transport cross section. The heuristic 
translation of these results for the sQGP (given in 48 ) 
requires the screening mass in the transport cross section 
to be replaced by 

Mscr — ► - 1-3T ~ gT. (17) 

Surprisingly this is -quantitatively- the same temperature 
dependence of the Debye mass as that obtained within our 
formalism (cf. Eq. (fTl"j) and discussion below). 

Some additional comments are in place: Whereas in 
our virial expansion approach the screening mass is calcu- 
lated from the thermodynamic quantities and the formal- 
ism includes crucial ingredients such as relativity and an 
effective screened potential - retaining thermodynamical 
consistency - the results of Ref. [JSJ are obtained within 
several approximations, i.e. a non-relativistic treatment, 
an unscreened electromagnetic interaction and (though 
motivated but) ad hoc modifications of the transport cross 
section. In this sense the latter approach is heuristic: the 
modifications do not automatically follow from the for- 
malism used but have been inserted ad hoc to achieve a 
better description of the experimentally findings. 



4 Results 
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Fig. 1. (Color online) Entropy density s of the QGP as a func- 
tion of the temperature divided by T 3 from the virial expansion 
(solid line) using g = 1.30. For comparison the corresponding 
SB limit is displayed by the dashed line. The 1QCD results 
(open squares) have been adopted from Ref. pQ. 

Before we discuss the results for rj/s from the virial 
expansion approach, we demonstrate the applicability of 
our formalism for the entropy density in comparison to 
the 1QCD calculations from Ref. [1 . In Fig Q] the entropy 
density s (divided by T 3 ) is shown as a function of the 
temperature expressed in units of the critical tempera- 
ture T c from the virial expansion approach using Eq. @ 
(solid line) as well as in the SB limit (dashed line). The 
symbols denote the 1QCD calculations from Ref. pQ. For 
completeness we show also the entropy density in the con- 
fined phase below T c , where we have calculated all ther- 
modynamic quantities (cf. Ref. [33]) within a generalized 
resonance-gas model and matched the different phases at 
equal pressure. Near T c the deviation of our results from 
the ideal gas limit are sizeable and huge in the confined 
phase. 

The ratio of the viscosity to the entropy density has 
been calculated using Eq. (|14l) . where the quark contri- 
bution of the entropy density s q is given by Eqs. ([T5]) 
and (|16|) . For the transport cross section the general ex- 
pression given in Eq. (jj) has been implemented, where 
the effective temperature dependent coupling ay(T) is 
given in Eq. (113[) and the Debye screening mass is given 
by Eq. (fTT]) . In Fig. [5] our main results are presented in 
comparison to other estimates. In the deconfined region, 
T/T c > 1, the solid red line show the results for r\j s as a 
function of the temperature (in units of the critical tem- 
perature T c ) using the EoS (|15H16[) . the coupling cvy (fT5|) 
and the Debye mass (jlip derived within the virial ex- 
pansion approach. Additionally, the experimental point 
(square) from [57] and the lattice data from Ref. [55] (tri- 
angles) and from Ref. [SB] (full dots) are shown for com- 
parison. In the confined phase, T/T c < 1, the purple region 
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Fig. 2. (Color online) The viscosity/entropy density ratio rj/s as a function of the temperature expressed in units of the critical 
temperature T c for T/T c < 1 and T/T c > 1. The lattice results are from Ref. [55] (triangles) and from Ref. [56 (full dots). The 
different lines denoted by xPT, RG and Virial exp. stand for the results from the scaling behavior of the chiral perturbation 
theory, the hadron resonance gas and the virial expansion approach (for T/T c > 1), respectively (see text). 



close to T c shows the estimates for rj/s in the resonance- 
gas model |58j , including all the known particles and reso- 
nances with masses to < 2 GeV and also an exponentially 
rising level density of Hagedorn states for to > 2 GeV. 
The dotted (red) line shows the scaling i)/s a T -4 com- 
bined with the requirement that rj/s = l/4ir at T c . This 
scaling behavior at low temperature has be found within 
chiral perturbation theory in Ref. [55] . Additionally, the 
range (0.8 — 1.5) for rj/s from perturbative QCD (pQCD) 
from Ref. [55] is sketched as a blue region. Furthermore, 
the KSS lowest bound is also indicated by the orange area. 

In the deconfined phase several features become evi- 
dent: 

i) At T c our result for rj/s ~ 0.097 is very close to the 
theoretical bound of l/(47r). Additionally, the convex- 
ity of the specific viscosity near T c suggests a minimum 
close to T c as expected in Refs. [381139] , 

ii) An almost linear increase of rj/s with the temperature 
is found for 1.5T C < T < 3T C . 

Hi) Currently the different lattice calculations are unable 
to provide quantitatively reliable results; the large er- 
ror bars of the lattice data do not allow for a conclu- 
sive comparison. Qualitatively, the increasing behavior 
of the specific viscosity with the temperature might 
be confirmed. In contrast, the experimental point (full 
square close to T c ) is reproduced very well by our result 
within the virial approach. 

iv) At higher temperatures a saturation of the ratio rj/s 
is found which roughly coincides within the band for 
pQCD. 

This last finding is supported by the observation that the 
entropy density s as well as the viscosity scale as oc T~ 3 



for high temperatures. For the entropy density s this is 
in line with the SB limit, which is approximately repro- 
duced by lattice pQ and model calculations [2Tfl[55] at high 
temperature. For the shear viscosity ij also very different 
approaches show this functional dependence: the strong 
quark-gluon plasma from AdS/CFT [BD], the quasiparti- 
cle approximation with differently modeled quark selfen- 
ergy pmit32"] as well as the weak coupling estimate from 
Ref. [H]. 

With respect to the confined phase the resonance-gas 
as well as the xPT calculation suggest a decreasing behav- 
ior of rj/s with increasing temperature. The 'constrained' 
scaling behavior (~ T~ 4 ) slightly underestimates the reson- 
ance-gas calculation from Ref. [58 . Note, however, that 
the entropy density in the Hagedorn resonance model sen- 
sitively depends on the level density employed for the con- 
tinuum states. Accordingly, we show these results in Fig. [5] 
by a shaded band. Choosing the scaling function (~ T~ 4 ) 
for the confined phase and our virial calculation for the 
deconfined phase we obtain a well connected description 
for rj/s in the whole temperature range (up to 5 T c ) which 
clearly shows a minimum near T c close to the KSS bound. 

To investigate in more detail the possibility of a min- 
imum of rj/s at the critical temperature from the virial 
expansion calculation we search for an approximation of 
the specific viscosity in the deconfined phase around T c . 
We employ a Taylor expansion of rj/s as a function of 
t = T/T c at t = 1, i.e. 

oo N 

^ = 1)" = 5> n (t-l)" + -Rj V) (18) 

n=0 n=0 
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Table 1. Results for the fit of r)/suptoT= 2T C 



N 


ai 


a.2 


03 


X 7dof 


2 


0.0968 


0.2359 




0.00102 


2 





0.3566 




0.03959 


3 


0.0616 


0.3528 


-0.08747 


5.1 10~ 5 


3 





0.5373 


-0.2163 


0.00255 



Table 2. Results for the fits of n/s up to T = 3T C . 



N 


0,2 


a3 ai 


xVdof 


2 


0.2477 




1.4478 


3 


0.4469 


-0.1194 


0.0164 


4 


0.5044 


-0.1998 0.0268 


0.0061 



with 



a " = ^!^F^ (< = 1) 



(19) 



In (fl~5|) N indicates the order of the Taylor polynomial 
while Rn stands for the rest of the corresponding expan- 
sion. Our strategy is now to investigate how different poly- 
nomial approximations - labeled by different N - can de- 
scribe n/s up to T = 2T C , T = 3T C and T = 4.5T C . For 
each N we calculate the x 2 /dof to evaluate the quality 
of the approximation. Since the expansion coefficients are 
connected to derivatives of the specific viscosity at T c (see 
Eq. (fF9")) ') by requiring that some coefficients vanish, we au- 
tomatically impose that the corresponding derivatives are 
equal to zero. In particular, with the requirement a\ = 0, 
we may investigate the possibility of a minimum at the 
critical temperature. Some polynomial fits to our results 
are given in the Tables [TJ [2] and [3l 

As expected, the power N = 2 cannot describe the 
behavior of n/s for T > 2T C . However, the inclusion of 
higher power allows to reproduce very well our data also 
for higher temperatures. For N — 3 and N = 4 we find 
a very good approximation for our calculated results for 
n/s up to T = 3T C and T = 4.5T C , respectively. Therefore, 
a power law approximation with a vanising or very small 
linear coefficient suggests a simple parametrization for the 
specific viscosity which includes a minimum close to T c . 



Ref. 33 . In this context the investigation of the interac- 
tion between the partons in the deconfined phase plays a 
crucial role to reproduce the thermodynamic properties of 
the QGP in comparison to lattice QCD. By using a gen- 
eralized classical virial expansion we have calculated the 
corrections to a single-particle partition function starting 
from an interaction potential whose parameters are fixed 
by thermodynamical quantities. Furthermore, in the virial 
expansion approach we can directly extract the coupling 
ay to be employed for the determination of the trans- 
port cross section which enters the ratio n/s. We find 
rj/s w 0.097 at T c which is very close to the theoretical 
lower bound. Furthermore, for T < 1.8T C the ratio is in 
the range of the experimental estimates 0.1 — 0.3 extracted 
from RHIC experiments. 

Additionally, a detailed analysis of the temperature 
dependence of our results for n/s has been performed. 
Within a Taylor expansion around the critical tempera- 
ture we found that a power law with a vanising or very 
small linear coefficient suggests a simple parametrization 
for the specific viscosity. This indicates the existence of a 
minimum in n/s close to T c . 

Since we focus on the deconfined phase only, we do not 
investigate whether in the vicinity of T c a phase transition 
or a rapid cross over occurs |63j . However, our approach, 
which provides a unified description of the QGP thermo- 
dynamic as well as of its sheer viscosity, is a first important 
improvement in the description of the quark-gluon plasma 
beyond mean- field models [55] . 

For further work, using relativistic molecular dynami- 
cal simulations, where the phenomenological quark-quark 
interaction used here can be implemented as well as the 
QCD equation of state, we may calculate further correla- 
tions in the partonic phase by considering partons in a box 
with periodic boundary conditions at fixed energy density 
(or temperature). Furthermore, such molecular dynami- 
cal calculations will allow to study partonic systems also 
out of equilibrium and provide important insight on the 
dynamics of hadronization. 
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Table 3. Results for the fit of n/s up to T = 4.5T C . 



N 


OL2 


03 


0,4 


xVdof 


2 


0.1227 






43.5251 


3 


0.3756 


-0.0845 




0.8524 


4 


0.4955 


-0.1781 


0.01738 


0.0228 



5 Conclusions 

We have performed an investigation of the specific vis- 
cosity n/s in the QGP in a dynamical way within kinetic 
theory using the virial expansion approach introduced in 
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